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Dynamical Chern-Simons gravity is an interesting extension of General Relativity, which finds 
its way in many different contexts, including string theory, cosmological settings and loop quantum 
gravity. In this theory, the gravitational field is coupled to a scalar field by a parity-violating 
term, which gives rise to characteristic signatures. Here we investigate how Chern-Simons gravity 
would affect the quasi-circular inspiralling of a small, stellar-mass object into a large non-rotating 
supermassive black hole, and the accompanying emission of gravitational and scalar waves. We 
find the relevant equations describing the perturbation induced by the small object, and we solve 
them through the use of Green's function techniques. Our results show that for a wide range of 
coupling parameters, the Chern-Simons coupling gives rise to an increase in total energy fiux, which 
translates into a fewer number of gravitational-wave cycles over a certain bandwidth. For space- 
based gravitational-wave detectors such as LISA, this effect can be used to constrain the coupling 
parameter effectively. 
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I. INTRODUCTION 

The inspirals of stellar- mass compact objects, typically 
black holes (BHs) or neutron stars, into supermassive 
BHs at the galactic centers are among the most promising 
sources for the space-based gravitational wave detector 
LISA 0. These processes arc known as Extreme Mass- 
Ratio Inspirals (EMRIs) and allow for stringent tests of 
general relativity to be done, for several reasons 0,0]. 
First of all, they emit 10^ cycles or more of gravitational 
radiation in the LISA band, during the timespan of the 
LISA mission (^1 — 5 years). Furthermore, this signal 
is emitted when the stellar mass object is close to the 
horizon of the supermassive BH, thus encoding the fea- 
tures of the strong-field BH spacetime and of the strong- 
curvature regime of general relativity Q • Finally, EMRIs 
are expected to be relatively clean systems, with negli- 
gible perturbations from surrounding matter. Therefore, 
the detection of the signal from an EMRI would allow 
us to test the strong-field regime of gravity, where possi- 
ble deviations from general relativity may show up. Of 
course, this would be possible only if we understand how 
alternative theories of gravity would affect the EMRI sig- 
nal. 
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In this paper, we discuss the imprint on EMRI sig- 
nals of dynamical Chcrn Simons (DCS) gravity, an ex- 
tension of general relativity [H-Q in which the Einstein- 
Hilbert action is modified by adding a parity-violating 
Chern-Simons (CS) term, that couples gravity to a scalar 
field. Among other proposed alternative theories, there 
are strong motivations to consider DCS gravity. Indeed, 
this correction arises in many versions of string theory 
and in loop quantum gravity [l0l - [l2j . Furthermore, 
CS gravity can be recovered as a truncation of low en- 
ergy effective string models [1, [Tsl a nd could also explain 
several problems in cosmology 

(see Ref. [13 for 

a recent review). 

DCS gravity would affect the EMRI signal in several 
ways. First of all, it modifies the spacetime metric of 
rotating BHs [l^l (though the metric of spherical, sta- 
tionary BHs is not affected). These modifications induce 
deviations in the motion of the stellar mass object, which 
affect the emitted gravitational signal. This effect has 
been studied in [2lj , where the deviation of EMRI orbits 
due to spacetime metric modification has been evaluated. 
In addition, in this theory the dynamical equations of 
gravity change, as the gravitational field is coupled with 
a scalar field. This is the effect studied in this paper, in 
which we evaluate the change of the EMRI signal due to 
the DCS modification of the dynamical equations of the 
gravitational field. 

We shall focus on the simplest case of EMRI system: 
the inspiralling of a test particle around a Schwarzschild 
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BH. Since the Schwarzschild metric is a solution of DCS 
gravity, the orbital effects studied in [2l| vanish. How- 
ever, as wc will show, DCS gravity could significantly 
affect the EMRI signal. We generalize the equations de- 
scribing perturbations of Schwarzschild DCS BHs, de- 
rived previously by some of us [22|,|23|, to include a source 
term describing an orbiting particle. Then, we solve these 
equations using a generalization of the Green's function 
method, evaluating the emitted flux of gravitational en- 
ergy and scalar energy. Finally, we determine how the 
modified energy loss affects the emitted gravitational sig- 
nal, by computing the change in the number of orbital 
cycles. Similar methods have been already successfully 
applied to stud y th e effects of neutron star structure in 
EMRI systems [2J], to investigate the nature of the cen- 
tral massive object [25j and also to put constraints on 
other alternative theories of gravity, like Brans-Dicke the- 
ory [H. 

We find that, although the DCS coupling only mildly 
changes the total energy fiux, the corrections accumu- 
late during the EMRI inspiral, producing a decrease 5M 
in the number of orbital cycles which is potentially de- 
tectable by LISA. Our results can be summarized by the 
following fit: 



- -26C\/— exp <j -1.21og^o 

m2 



mi 



(1.1) 



where mmax ^ IO^Mq (its precise definition is given in 
Eq. (14.21^ ). C is a parameter characterizing the DCS cor- 
rection, and mi, 7712 are the masses of the supermassive 
BH and of the stellar mass object, respectively. The fit 
above is valid in a wide region of the parameter space, 
i.e. mi e [10^ 10^]Mq, m2 S [1, 10]Mo and C < 1, with 
an error of at most a few percent. For larger values of C, 
the fit (|l.ip still gives an order-of-magnitude estimate of 
6N: for 1 < C < 20, the error does not exceed 50%. 

The plan of the paper is the following. In Section [H] 
we discuss the equations describing perturbations of the 
Schwarzschild metric in DCS gravity induced by a point- 
like particle on a circular orbit about the BH. In Sec- 
tion UTT] we describe the Green's function approach to 
solve the perturbation equations. In Section HV Al we de- 
rive the energy flux associated to the gravitational and 
scalar radiation, and compute this fiux solving numer- 
ically the perturbation equations. In Section IIVBI we 
determine how this energy flux affects the gravitational 
signal. In Section|V]we draw our conclusions. The deriva- 
tion of the perturbation equations with source is dis- 
cussed in detafl in Appendix 1X1 Appendix IB] is devoted to 
describe a perturbative Green's function approach, which 
is valid in the small coupling limit and may be potentially 
useful for future analytical calculations. In Appendix [C] 
we compare different prescriptions to compute the num- 
ber of orbital cycles. 



II. DCS GRAVITY AND PERTURBATION 
EQUATIONS 

The action of DCS gravity reads [1^ (we use geomet- 
rical units c = G = 1) 

S^ — I d^x^R + - [ d'^xyf^'d*RR 
16tt J ■ A J 

-^J d^xV^ [g^'Wa^^b^ + V{{))] + 5,„at , (2.1) 
where i? is the scalar field and 

* T>r> n * Tjbacd ^ t3 ^baef ncd in r)\ 

^^—n.abcd n ~ -^^abcd^ ef ■ l^-^j 

Neglecting the scalar potential V{'d), the equations of 
motion are 

Rab = -IGTraCab + (Xab - ^gabTj (2.3) 
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'RR 



(2.4) 



where the stress-energy tensor accounts for the matter 
and the scalar field contributions. Tab = T^b'^^ + T^f^, with 



and 



Tab = ( ^■,a^;b - ^^:.c^--' 



fjab ^ ^.^eCde(ay^^fc)^ _^ ^^^^ * j^d{ab)c 



(2.5) 



(2.6) 



Since in any spherically symmetric background * RR = 
and C"'' = 0, spherically symmetric solutions of general 
relativity, like the Schwarzschild metric 



ds^ = -fdr + f-'dr^ + r^dO.^ , 



(2.7) 



arc solutions of DCS gravity [19|. Here, M is the ADM 
mass of the spacctime and 



/ = l-2Af/r. 



(2.8) 



A remarkable feature of EMRIs is that they can be 
described with great accuracy within a perturbative ap- 
proach and in an adiabatic approximation. Indeed, dur- 
ing most of the inspiral the stellar-mass object can be 
considered as a test particle moving in a single massive 
BH background, the timescale for merger being much 
longer than a single orbital period. Hence, at each in- 
stant, we consider that the particle follows a geodesic of 
the BH spacctime and the geodesic parameters, i.e. the 
orbital energy and angular momentum of the particle, 
would change adiabatically: they can be computed by 
solving the linearized Einstein's equations for geodesic 
motion. In this way one finds the inspiralling orbit and 
the corresponding gravitational waveform, as explained 
below. This procedure takes into account the main effect 
of the back-reaction (the so-called "non-conservative part 
of the self-force" ) . A more detailed analysis, which would 
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also consider the "conservative part of the self- force" (see 
m HI for the case of non-rotating BHs and [H, for 
a review), is beyond the scope of this work. 

Remarkably, due to Eq. (|2.4p . the tensor (|2.6|) satis- 
fies yCah = SttVT^j. Using the latter equation and 
the Bianchi identities, the matter stress-energy tensor is 
conserved also in DCS gravity, i.e. VT^''* = 0. It fol- 
lows that point-like particles travel on geodesies, exactly 
as in general relativity. Furthermore, we assume that 
the point-particle is non-spinning, neglecting the effects 
of spin-orbit interactions. These effects are proportional 
to the mass-ratio and they can be safely neglected in the 
study of EMRI systems [3l|. 

We shall study the perturbations of a static, spherically 
symmetric BH of mass M, due to a non-spinning point- 
like particle of mass ^ on a circular orbit around the BH. 
As discussed in [2^, the Schwarzschild metric (|2.7p with 
a vanishing scalar field is the only static, spherically sym- 
metric BH solution in DCS gravity. In this background, 
we expand the gravitational and scalar perturbations in- 
duced by a point-like particle in tensor spherical har- 
monics, building the Zerilli and Regge- Wheeler functions, 
Z^"(r), Q^™(r), and the scalar field function e^'"(r). In 
the frequency domain, the perturbation equations read 
(see Appendix [Xl for details): 



dri 



dr'l 
dri 



0^"(r)=rw(r)e^™(r)-l-5lJ|V(r) 

(2.9) 

e^'"(r)=T5(r)g^'"(r) + 4"(r) 



(2.10) 
(2.11) 



Z'^^ir) = S'|'"(r). 



In the equations above r* is the tortoise coordinate de- 
fined by dr/dr^ = f, and the potentials read 



1) 6M \ 



TRwir) = f 



96i'KMuja 



Vsir) = f 
Tsir) 

Vz{r) 



\ r2 



B/3 



2M 



-I 
f 



{e + 2y.6Mia 

{e~2y. r^puj 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



2A^ A+1- 



r 



IMP 



A+ — 
r 



(2.16) 



where A = (£ + 2)(^ - l)/2 and A = A + 3M/r. The 
source terms are given in Appendix El and (cf. 
Eqs. (jA10[) and (IA14 [ )) are the same as those computed 
in general relativity [32|, whereas S'f"" (cf. Eq. (|AT21) ) is 
proportional to a and it is a novel term introduced by 
the CS coupling. 



We stress that Eqs. p.9|) and p.lOp are coupled 
through the CS coupling a. In the sourceless case (i.e., 
no exterior matter), they reduce to those in Refs. p^[23| 
with Q^™ = zw^"^™ (notice the sign difference due to 
our definition (|A8[) ). In the general relativity limit 
a = 0, they decouple into the Regge- Wheeler equa- 
tion with sources [32 . [33| plus Klein-Gordon equation 
without source. Interestingly, in the DCS Schwarzschild 
background the polar (even parity) gravitational sector, 
Eq. (|2.1ip . decouples from the scalar sector. Thus the 
CS coupling does not affect the Zerilli equation, which 
simply reads as in general relativity js^ . 

To conclude this Section, we remark that if we rescale 
the scalar field, in order to express its kinetic term in a 
canonical form 



(2.17) 



then the DCS action takes the form 

S = — f d^x.r^R + - ^ 
167r 

1 

^2 



16V^ 



d^x 



(2.18) 



where, following Ref. [20], we have defined ^ = IQira^/fi. 
Therefore, we expect that all physical observables depend 
on the parameter ^ (see also the discussion in [1^), or 
equivalently on the dimensionless parameter 



c 



167ra^ 

/3AP 



(2.19) 



where M is a quantity with the dimensions of mass (i.e. 
of length) associated to the system under consideration 
(in our case, the mass of the supermassive BH). This is 
indeed the case, as we shall show, for the motion of a 
test particle around a static, spherically symmetric BH. 
This procedure, of rescaling the scalar field to have a 
canonical kinetic term in order to get rid of redundant 
parameters, is well known in the context of scalar-tensor 
theories (see [s^] and references therein). However, here 
we shall follow the formulation of Ref. [201 j where a and 
/3 are kept as independent parameters. 



III. GREEN'S FUNCTION APPROACH 

In order to compute the gravitational-wave emission 
of a particle in geodesic motion around a spherically 
symmetric BH in DCS gravity, we shall solve the equa- 
tions (|2.9p - (|2.1ip by extending the standard Green's 
function techniques. As explained below, we work out the 
basic equations for general orbits, and then we specialize 
to circular motion. Our results can be easily generalized 
to eccentric orbits. In Appendix [B] we develop a per- 
turbative approach, which is valid in the small-coupling 
limit, and we compare it with the general method dis- 
cussed in this section (the agreement is very good for 
small coupling). 
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Even sector 



We start by considering the Zerilli equation (|2.1ip . 
which is not modified in DCS gravity. One considers two 
solutions Z^"- of the associate homogeneous equation 



2 



Z+ = 0. 



(3.1) 



(hereafter, we leave implicit the £, m indices) such that 



Z4 



e 



±0O . 



(3.2) 



Then the general solution reads 



Z(r) 



1 



Z+{r) / dr,Z-Sz + Z-{r) 



+ 00 



dr^ Sz 
(3.3)" 

where Wz = f{Z-Z'_^ - Z+Z'_) is the Wronskian 
and the prime denotes derivative with respect to the 
Schwarzschild radial coordinate, r. At infinity and at 
the horizon, where the energy flux is computed in term 
of Z, we get 



Z{r^ 



±oo) 



Wz 



dr^ZzrS 



(3.4) 



As we show in Appendix \X[ for a circular orbit at r = f 
the source term has the form Sz ^ S{r ^ f). In this case 
to compute the integral p.4|) it is sufficient to evaluate 
the integrand at r = f. 



B. Odd sector 



Let us consider the modified Regge- Wheeler equa- 
tion (|2.9p . coupled with the scalar equation (|2.10p : 



^ + - - Vnw 



Q = Srw H ^ aQ , (3.5) 



Q = Ss-f 



{£ + 2)l6iMa 



(3.6) 



Introducing P^" = dQ'^"'/dr^ and = de'^'^'/dr^, we 
can write the equations above as a first order system 



= S. 



(3.7) 



where * = (g,e,P,$)^ and S = {0,0, Srw, SsV are 
four dimensional vectors and V reads 



V 









-1 











- Vrw 


-Trw 





Ts 








(3.8) 



The system p.7p can be solved by standard methods (see 
e.g. Ref. d^). For this purpose, define the 4x4 matrix 



X whose nth column contains the nth solution of the 
homogeneous system dx/dr^, + Fx = 0, i.e. Xij = 
where the j index denotes a solution of the homogeneous 
system and i is the vector index. It can be shown that 
also the matrix X constructed in such a way is a solution 
of the associated homogeneous system, in the sense that 



dX 

dr^, 



VX = 0. 



(3.9) 



In order to solve p.7p . we impose the ansatz ^ = Xa, 
where H is a vector to be determined. Substituting the 
equation above into the inhomogeneous system and using 
Eq. dSH) we find 



X'^S, 



(3.10) 



(3.11) 



The matrix X, contains four independent solutions of 
the homogeneous system, supplied by suitable boundary 
conditions. We impose 



and the solution to p.7p reads 

* = X / d7\X-^S. 



A± 
B± 



zbcX) . 



(3.12) 



As explained in Ref. |2j], two linear independent solu- 
tions can be constructed by choosing: (i) A± = 1 and 
an arbitrary B± = B^^^ and (ii) B± = 1 and an arbi- 
trary A± — A^£\ provided A^^^ B^f' ^ 1. This procedure 
can be applied twice: first we construct two solutions, 
, imposing boundary conditions at the hori- 
zon and integrating outward, and secondly we construct 
two further independent solutions, Ix^^x^'l, by im- 
posing boundary conditions at infinity and integrating 
backward. 

Finally, from Eq. p. lip , we can write the solutions for 
the gravitational and scalar waveform, which satisfy the 
correct boundary conditions, as follows 

2 
i=l 

©W = E {(^fir)I^\r) + e«(r)/«(r)) ,(3.14) 

i=l 

where 

/i^^ = £ + D^I>Ss) , = 2) 

and the functions Cj.'"* and Dj!'' depend on the so- 
lutions of the homogeneous system, ^xL^^x^^l and 

Ix^i^' , x^^^ |, and can be straightforwardly computed 

from the components of the vector X~^S. For complete- 
ness, their expressions are given in Appendix [P] 
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Circular orbits 



IV. RESULTS 



So far our approach generically holds for a point-hke 
particle in geodesic motion. However, the above formulae 
simplify significantly in the case of circular geodesies at 
r = f, as we show in Appendix |^ In particular, the 
source term for the scalar equation vanishes, Sg = 0. 
The remaining source terms can be factorizcd in order 
to extract a Dirac delta contribution ~ i5(a; — niiOK), 
where luk is the Keplerian frequency (jA26p . In the rest 
of this Section we write explicitly the indices £, m and 
the dependence on uj. 

To compute the Zcrilli function we replace Eq. (|A33|) in 
Eq. (|3.3p . Integrating by parts to get rid of the derivative 
of the delta function, we get Z^"'(a;,r) = Z.^™ 
muix) with 



(r) 



Wz 



/ 



for r ^ f respectively; hatted quantities arc defined in 
Appendix |X1 to which we refer for further details. At the 
boundaries, if wc call ZjP^{uj) = Zj!"(tj,r* ±00), we 
obtain 



f \ f 
For later use, we will write this as 



S{uj—mujK)e 



(3.15) 



where is a constant. 

A similar procedure can be applied to the axial sector. 
Using Eq. (jA32p . the Fourier transform of the Regge- 
Wheeler and scalar function, at the boundary r» — >■ ±00, 
read 



em , 



05" (cj) = 0':^s{uj 



(3.16) 
(3.17) 



where, from Eqs. (I3.13P and p.l4[) . we have 




As we show in the next section, the energy flux at infin- 
ity and at the horizon can be computed in terms of the 
quantities Z^™, Q^™ and 0^™ given in Eqs. (|XT51) . 
and (|3.17p . respectively. 



A. Energy flux 

The flux of gravitational energy can be computed in 
terms of metric perturbations. At asymptotically flat, 
future, null infinity, the expressions of the effective grav- 
itational wave stress-energy tensor (i.e. the Isaacson ten- 
sor) in DCS gravity and in general relativity coincide 
[36| . Therefore, we can use the machinery derived in the 
framework of general relativity to determine the emitted 
gravitational energy fiux at infinity in DCS gravity. 

On the other hand, the expression of the rate of energy 
absorbed by the horizon is also formally equivalent to 
that in general relativity (sec e.g. Ref. [331). Indeed, the 
derivation involves the first law of BH thermodynamics, 
which relates the change in energy M with the change 
in the horizon area A, and Raychaudhuri's equation to 
calculate A in terms of the shear tensor. In DCS gravity, 
both these steps proceed exactly as in general relativity, 
all the dynamical information being eventually encoded 
in the waveforms. 

Hence, the energy fiuxes at (null) infinity and at the 
horizon formally read as in general relativity [ssj 



dEr, 



dx 



647r (£ - 



2)! ^ 



|ii"(x)p-f4|Qi"(x)p 



(4.1) 



where the sum is taken over negative and positive m and 
X = t^r* are the retarded and advanced coordinates, re- 
spectively. The inverse Fourier transform of p.isp . p.l6p 
is 



dLLiZ±S{uj — mwK)^ 



dujQ±5{uj - ma;K)e-'""^-(*T'^*) 



therefore 

E: 



(4.2) 

On the other hand, the scalar energy fiux reads (see, 
e.g. Ref. IH) 



Escai = ~r^fir) / dnr^ 



seal 



(4.3) 



From the stress-energy tensor of the scalar field, T^l^^ = 
/3(V(ai?*V,,)i?-l/25Qf,Vci?V'=i?*). Inserting Eq. ^M^ and 
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using the asymptotic behavior at infinity (|3.17p . the en- 
ergy flux reads 



p± = 

seal 



dE. 



seal 



dx 



Im 



(4.4) 



Finally, since the orbital frequency is related to the 
orbital velocity v and to the semi-latus rectum (which 
for circular orbits is simply p = f/M) hy the relations 



iMLJKr=p 



1/3 ^ „-l/2 



(4.5) 



the energy flux E can also be considered either as a func- 
tion of V or p. The condition for the existence of stable 
circular orbits, f > nsco — 6M, constrains the values of 
V and p to p> 6 and v < 6~^/^ ~ 0.408. 

The method described above has been implemented 
in Mathematica. In our numerical approach, we have 
considered a series expansion at the horizon and at in- 
finity up to order eight for the boundary conditions p.2p 
and p.l2p . Our results arc summarized in Fig.[TJ When 
C = 0, our results agree with those in general relativ- 
ity [H,!!^] within one part in 10^ or better. Furthermore, 
in the small ^ limit, we develop an independent method 
(discussed in Appendix [B]) whose results are in perfect 
agreement with the ones discussed here. 

As expected, the CS corrections are more effective 
when p ~ 6, i.e. close to the innermost stable circu- 
lar orbit (ISCO), where circular orbits probe the strong 
curvature region around the massive BH. Far away from 
the source the CS contributions are negligible. This is 
clear from the left panel of Fig. [TJ where we show the 
relative difference in the emitted power 



SE 



E, 



DCS 



-E, 



GR 



E, 



GR 



E, 



GR 



where Edcs 



E" + E°° - 

grav ' grav 



H 



rpoo 
' ^scal ' 



(4.6) 



i.e. it is the 



sum of the contributions coming from the gravitational 
and scalar fluxes, both at the infinity and the horizon, 
and Eq^ = ii'DCs(C = 0), i.e. the energy fiux in general 
relativity. Clearly, the scalar contribution to -Egr is van- 
ishing. The relative difference is positive, i.e. the total 
power emitted in DCS gravity is larger than in general 
relativity. This is consistent with the fact that in this 
theory there is an extra scalar degree of freedom, which 
introduces further energy dissipation channels. Although 
the difference in the total flux is positive, we find that for 
some subdominant {£ > 3) mode, the energy flux may be 
smaller than the corresponding flux in general relativity. 
This shows that a conversion of scalar into gravitational 
energy is possible, due to the CS coupling. 

Furthermore, even if the axial flux can be as large as 
twice the axial flux in general relativity (for example 
when p ~ 6 and for ( ~ 10), the correction to the to- 
tal energy flux is significantly smaller. Indeed, the lead- 
ing contribution to the energy fiux arises from £ = m 
modes, which have polar parity because selection rules 



imply that for even values of £ + m only polar perturba- 
tions are sourced (see Appendix [X\ for details). There- 
fore, since DCS corrections only affect the axial sector of 
Schwarzschild perturbations, their contribution is sub- 
leading with respect to that coming from the polar per- 
turbations. In the most favorable case (p ^ 6 and ^ ^ 10) 
the total energy flux (summing over polar and axial con- 
tributions up to £ = 5 and ^l < m < I) only differs from 
the general relativity value by a few percent. Typically, 
the deviation is smaller, as shown in the left panel of 

Fig.m 

In the right panel of Fig. [1] we show the four contribu- 
tions to the total emitted power. Remarkably, the main 
contributions arise from the gravitational and scalar flux 
at the horizon, which are positive and sensibly larger than 
the contributions at infinity; see Appendix |B] for a dis- 
cussion on this behaviour. 

Note that the correction to the gravitational flux at 
infinity is negative but, since this is a subleading contri- 
bution, the correction to the total energy flux is never- 
theless positive. 

Finally, from the results shown in Fig. [T] we can extract 
the following dependence in the small v limit 



SE, 



7H 

grav 



SE^ 



,12 



-^GR 



,14 



(4.7) 
(4.8) 



Notice that the high power in v introduces large errors in 
the fits above and these results should be understood as 
lower limits for the post-Newtonian (FN) orders of the 
CS effects. The scalar flux at infi nity is consistent with 
analytical predictions at FN level |41| . 



B. Gravitational- wave signal 

Although the difference between the energy emission 
by EMRIs in general relativity and their emission in 
DCS gravity is small, stellar-mass objects can orbit su- 
permassive BHs for ~ 10^ cycles (many of which can 
occur near the ISCO) while in the sensitivity window of 
LISA, before reaching the ISCO and eventually plunge. 
Hence, the small deviations accumulate and they can re- 
sult in sensible modifications when one looks at the entire 
inspiral-t-merger process. Let us give a rough estimate of 
this effect, and of its observability by space based detec- 
tors like LISA. 

A useful quantity to consider is the number of 
gravitational-wave cycles accumulated within a certain 
frequency band (see for instance Rcf. 0|). This quan- 
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FIG. 1. (Color online) Left: Relative difference between the power emitted in gravitational waves in general relativity and DCS 
gravity, cf. Eq. (|4.6p . for different values of (. The sum is truncated at ^ = 3. The contribution in DCS gravity includes both 
the gravitational and the scalar fluxes, -Edcs = E^rav + iSscai at infinity and at the horizon. Right: Different contributions to 
the relative difference for different values of ^ = 0.5, 0.4, 0.3, corresponding to different colors. The main contributions arise 
from the fluxes at the horizon and they are positive, whereas the difference in the gravitational flux at inflnity is negative. 



tity is defined 

Af= f^' Idf, (4.9) 
Jf, J 

with 

fi = max (/low, Ayr) , (4.10) 

// =mm(/isco,/up) • (4.11) 

Here /low = 10"^ Hz and /up = 1 Hz arc two typi- 
cal (lower and upper, respectively) cutoffs of tlic LISA 
noise curve, /isco = (6'^/^7rM)~^ is the frequency at the 
ISCO. Finally, /lyr is the frequency of the gravitational 
radiation emitted one year before the coalescence (coales- 
cence is for simplicity defined to occur when the particle 
reaches the ISCO). This choice is due to the fact that, 
despite the timespan of the LISA mission will likely be 
~ 1—5 years, we make the most conservative assumption, 
considering one year of observation time. The frequency 
/lyr is obtained by solving 

(•/isco 

Tobs = / ^ = 1 year. (4.12) 

"'/lyr / 

To solve Eqs. (|4.9p and (|4.12l) . one has to know how the 
particle is inspiralling into the massive BH. We use a sim- 
ple prescription, the so-called adiabatic approximation. 



in which the back-reaction is obtained by the energy flux, 
that is computed assuming that the particle moves on a 
geodesic orbit. Later in this Section we shall discuss this 
approach in more detail, assessing its accuracy; here we 
only remark that the error due to this approximation is 
much smaller than the other uncertainties in our problem 
(specifically the observation time, but also fundamental 
issues such as the magnitude of the coupling constants). 

In order to compute /lyr using Eq. ()4.12|) . one may con- 
sider a PN formula for /. This choice has been widely 
adopted in previous works [H, H^,!!!]- For EMRIs, the 
velocity during the latest stages before coalescence can 
be as large as a fraction of the speed of light and the PN 
approach may introduce sensible errors, as discussed in 
Appendix [Cl Therefore, we adopt a different prescrip- 
tion, as explained below. Our approach is well suited for 
the large-velocity and strong-field regime and thus it is 
expected to be more accurate for EMRIs. In Appendix 
Owe compare it with other, less accurate, approaches. 

We consider geodesic motion of a two-body system 
with masses rrii and m2. The motion can be effec- 
tively reduced to that of a particle with reduced mass 
fi = mim.2/M ~ rjM orbiting a central object with mass 
M = TOi +m2, consistently with the notation used in the 
previous sections. The frequency of gravitational waves 
emitted by the particle on a circular geodesic at r = f 
reads 



^ In this section / indicates the gravitational wave frequency, 
whereas, in the previous sections, it denoted the {t, t} compo- 
nent of the metric 1 12. 7| I. To avoid confusion while keeping the 
standard notation, in Eq. 114.1311 this component is denoted by 
900(f)- 



UJK 



from which we get / 



2f 



(4.13) 



\f ^. For a particle in a circular 
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TABLE I. Corrections to the number of gravitational wave cycles accumulated within the frequency band / £ [fi, //] for some 
typical two-body systems. Data correspond to Tobs = lyr and are roughly fitted by Eq. (|4.20[l . 



M 


(1.4+ 1O'')M0 


(1.4 + 7 X 10*)Mq 


(1.4 + 4 X 1O'')M0 


(10 + 4 X 10")Mo 


(10 + 1O'')M0 




Ji(nZ) 




n 0205 




0.0041 




0.0010 




00088 




00040 




//(Hz) 




0.4396 




0.0063 




0.0011 




0.00110 




0.00044 






—SN' 


—5M/N 


~5M 




—SM 




-5M 




—57V 


-SM/M 






38.2 


3 78 X 1 fl~^ 

O. (OA. LXJ 


900 


6 19x1 0~'^ 


455 


1.40 X 10~^ 


295 


q 94 X 1 0~'^ 


179 


1 38 X 10" 

± . <JO /N i-\J 


2 


1 n 

lU 


15.2 


1.51 X 10^^ 


365 


9 48 X 1 0"'^ 


188 


78 X 10"^ 


121 


4.07 X 10"'^ 


73.9 


5 69 X 1 0" 


3 


o 




R 91 X 1 fl^" 


152 


1.04 X 10"'^ 




2.44 X 10"^ 


iJU.O 


1.71 X 10"'^ 


31.2 


Zi.'-±\J A _LU 


3 


9 


2.12 


2.09 X 10"" 


51.7 


3.51 X 10"* 


27.2 


8.34 X 10"* 


17.3 


5.82 X 10"* 


10.7 


8.20 X 10" 


4 


1 


0.98 


9.70 X lO"'^ 


24.0 


1.64 X 10""* 


12.7 


3.89 X 10"* 


8.05 


2.71 X 10"* 


4.97 


3.82 X 10" 


4 


0.5 


0.47 


4.65 X 10~'^ 


11.5 


7.85 X 10"^ 


6.09 


1.87 X 10"* 


3.87 


1.30 X 10"* 


2.39 


1.84 X 10" 


4 


0.4 


0.37 


3.69 X 10-'' 


9.15 


6.23 X lO"'' 


4.83 


1.48 X 10"* 


3.07 


1.00 X 10"* 


1.89 


1.45 X 10" 


4 


0.3 


0.28 


2.74 X lO"'^ 


6.81 


4.63 X 10"^ 


3.59 


1.10 X 10"* 


2.28 


7.68 X lO"'' 


1.40 


1.08 X 10" 


4 


0.2 


0.18 


1.81 X 10"'^ 


4.50 


3.06 X 10"^ 


2.37 


7.29 X 10"* 


1.51 


5.08 X 10"^ 


0.93 


7.17 X 10" 


5 


O.f 


0.09 


8.97 X 10~* 


2.23 


1.52 X 10"^ 


1.18 


3.61 X 10"^ 


0.75 


2.52 X lO"'^ 


0.46 


3.55 X 10" 


5 


0.01 


0.009 


8.95 X 10"^ 


0.221 


1.51 X 10"'' 


0.12 


3.61 X 10"" 


0.074 


2.50 X 10"" 


0.046 


3.53 X 10" 
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FIG. 2. (Color online) Left: Total number of cycles A/" as a function of the central mass mi, for several values of 7712, Tobs = lyr 
and C, — 0.5. Right: Corrections to the number of gravitational wave cycles accumulated during the inspiral of a small object 
(m2 — 1.4 Mq) around a supermassive BH of mass mi in 1 year observation time before coalescence. We show as a function 
of mi for different values of C^. Same legend as in Fig. [1] 



orbit we have 



f-2M 



^f{f - 3M) 
Putting all together, wc find 

/_ 3/ df ■ 

2 r dhnrh 



Finally, using Eq. (|4.13p wc get 

3 [Ml/3/-2/3_ 3^^2/3] 3/2 



(4.14) 



(4.15) 



7r2/3^^M5/6 r^,^i/3y-2/3 _ 6M7r2/3l 



f2/3 



DCS 



(4.16) 



where E'dcs = ^-grav + ^'scai = —Eorb IS the total en- 
ergy flux radiated away, using the standard flux balance 
equation 



En 



E. 



0. 



(4.17) 



In the balance equation above, both the contributions at 
the horizon and at inflnity must be taken into account. 

We remark that we are modeling the EMRI orbit in 
the adiabatic approximation: the particle is in nearly 
geodesic motion, allowing to compute, at each time, the 
emitted energy flux Ejjcs assuming a geodesic orbit; fur- 
thermore, in this approach one assumes that the back- 
reaction is described by the flux balance (|4.17p . This 
approximation neglects the so-called "conservative part 
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FIG. 3. (Color online) Correction to the number of gravita- 
tional wave cycles accumulated during the inspiral of a small 
object around a supermassive BH of mass mi — 4 x 10^ 
in 1 year observation time before coalescence, as a function 
of the mass m2 of the small object and different values of 
(same legend as in Fig. [1]). 



of the self-force" , which is a higher order effect contribut- 
ing marginally to the gravitational signal (see Sec. 2 of 
j44 | ) . Indeed, as shown in [l^ , for spinning BHs the con- 
servative contribution can account for at most few cycles 
of the entire process, while for non-spinning BHs it con- 
tributes less than one radiant to a one year evolution [i^ . 
Therefore, it may be necessary to take it into account in 
the data analysis of the process, but this effect can be 
neglected in assessing the relevance of DCS corrections 
for the EMRI signal. 

Using the energy fluxes computed in the previous Sec- 
tion, we can therefore obtain the frequency at any in- 
stant prior to merger and the number of cycles Af left to 
merger. The total number of cycles is shown in the left 
panel of Fig. O 

Wc can also compute the correction to the total num- 
ber of cycles due to the CS coupling as follows. From 
the relation (|4.16p . we obtain Sf/f = SE/Egr., with 
SX = Xdcs — -^GR- Then, at first order. 



/.J 
U f 



1 



SE 
Egr 



df, 



and the correction in the number of cycles reads 

rff f SE 

>f' f Egr 



SAf_ 



(4.18) 



(4.19) 



The difference in the number of cycles is shown in Table U 
for different values of the CS coupling ( and for a typical 
set of parameters. 

The corrections depend on mi, TO2 and on the CS 
coupling This is shown in the right panel of Fig. [5] 



and in Fig. |3l where we indicate the fiducial threshold 
SAf = {2tt)~^ cycles (i.e. (5$ = 1 rad) with a horizon- 
tal line. Corrections to general relativity are generally 
considered significant if they exceed one radiant over the 
observation time [13] . 

In the right panel of Fig. [3] we show the dependence 
on 7712 and on C for a central supermassive object with 
nil = 4 X 10^ Mq (i.e., the mass of the supermassive 
BH at the center of the Milky Way (Hi). Note that the 
dependence on m2 appears to be very mild. 

Overall, our results are well described by 



6J\fr 

where 



-26C 



Me 

TO2 



exp < —1.2 log 



10 



6.6 X lO^Mr. 



TO2 



(4.20) 



(4.21) 



is the location of the maximum in the left panel of Fig. [5] 
and it does not depend on the CS coupling. The fit above 
has been inspired by the curves in the left panel of Fig. [5] 
In a semi-logarithmic scale, these curves are approxi- 
mately Gaussian, y = yo exp [ao{x — Xm)^] , where yn, ao 
and Xm are the fit parameters and the expression (|4.20p 
is simply written in the coordinate x = logiQ(mi). We 
estimate an error on the fit smaller than a few percent 
when mi e [10^, 10^]Afo, m2 € [1, 1O]M0 and C < 1- 
As shown in Fig. |31 for larger values of the CS coupling, 
\SAf\ grows faster than linearly as a function of C, and 
Eq. (|4.20[) would acquire higher order in contributions. 
For example when C ~ 20, the fit (|420)) is accurate within 
50%. 

The presence of a maximum in SM{mi) (Fig. [21 right 
panel) can be understood as follows: on one hand, the 
DCS correction for a given value of C = Wna^/ipM'^) 
becomes more significant as M = nii + m2 increases, as 
it appears from Eqs. (|2.9p - (|2.16p : on the other hand, for 
large values of M the total number of cycles decreases 
(see Fig. [2l left panel), and thus SAf decreases, too. 



V. CONCLUSIONS 

We have studied the gravitational wave emission by a 
small object on a quasi-circular geodesic around a static, 
spherically symmetric, massive BH, in the context of Dy- 
namical Chern-Simons gravity. This process can describe 
for instance the inspiralling of a neutron star or stellar- 
mass BH into a supermassive BH, and is thought to occur 
frequently in the universe. In fact, EMRIs are one of the 
main preferred sources of gravitational waves for space- 
based detector LISA. We have shown that, because the 
stellar-size object spends many cycles in the bandwidth 
of LISA, the small effect of coupling to the DCS term 
"piles up" giving rise to measurable effects, in particular 
a decrease in the number of cycles over a fixed frequency 
bandwidth. 
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FIG. 4. (Color online) Correction to the number of gravitational wave cycles accumulated during the inspiral of a small object 
around a supermassive BH of mass mi in 1 year observation time before coalescence. Left: as a function of C, for some 
values of mi and m2 = IAMq. Right: same with m2 = IOMq. The corrections are linear in C,. 



Extensions of this work are necessary before useful con- 
straints to the theory can be obtained. In particular, the 
space of Schwarzschild BHs is a set of measure zero in the 
space of solutions: real BHs are most likely to be rotat- 
ing. An extension of the present formalism and results to 
spin ning; BHs (for instance to the general class found in 
Ref. [401) is highly desirable, but might have to wait for 
more powerful techniques or full-blown numerical simu- 
lations. Since the curvature invariants are higher close 
to rapidly spinning BHs and since the ISCO gets closer 
to the horizon, rapidly spinning BHs are a potentially 
very interesting tool to test DCS gravity. Furthermore, 
we have focused on circular orbits, but eccentricity may 
play an important role. The formalism we developed al- 
ready allows for studies of eccentric EMRIs in DCS grav- 
ity. For eccentric orbits around rotating BHs much larger 
corrections due to the CS coupling are expected. There- 
fore, our results should be seen as a lower limit of DCS 
corrections to the EMRI gravitational-wave signal. 

Finally, the results described in this paper are valid 
for any value of However, astrophysical observations 
already constrain the CS coupling [201 



(5.1) 



Notice that a and j3 are the physical parameters entering 
the action (|2.1|) , so that the constraints on ^ (for instance 
those possibly arising from a comparison between future 
observational data from LISA with Table U or Figs.[2]|4]) 
should always be converted into constraints on ^, depend- 
ing on the lengthscale of the system under consideration. 

In this case, from the definition (|2.19p and the con- 
straint (|5.ip we obtain that, for objects with mass M < 
IO'^Mq, the CS coupling ( may be larger than unit and 
our results show that in this case the CS coupling intro- 
duces large corrections, which are potentially detectable 



and not already ruled out by previous constraints. How- 
ever, the point-like approximation used in this paper is 
justified for EMRIs (mi/m2 ^ lO'' or larger) whereas, 
for smaller values of the mass ratio the structure of the 
smaller object, and the details of its backreaction, should 
be taken into account. 
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Appendix A: Derivation of the perturbation 
equations 

In this appendix we describe the derivation of the per- 
turbation equations (|2.9p . (|2.10p and (|2.1ip with source 
terms Ss and Sz respectively. 

Let us start with some definitions. We denote the 
Regge- Wheeler (RW) function by Q and the Zerilli func- 
tion by Z. F^™ are the usual scalar spherical harmonics, 
in term of which, following Refs. [s^lssj we define vector 
and tensor spherical harmonics as follows 



X 



B 1 



U' 



AB 



UabY 



'^AB — ^1 



\AB 



UabY 



w 



AB 



^(A\B) 
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where flAB = (l,sin^0), a bar denotes the covariant 
derivative with respect to the metric Qab , and Eab is the 
Levi-Civita tensor on the unit two-sphere. Hereafter cap- 
ital roman indices run over the angular coordinates (9, tp), 
while lower-case roman indices run over t and r. We 
work in the frequency domain and all the quantities are 
intended as Fourier transforms of some time-dependent 
quantity, i.e. schematically A{uj^r) = ^ J dtA{t,r)e^'^* . 
Finally we define A' = dA{uj,r)/dr. 



odd parity perturbations read 

lujh' + h'' + hi i /lo 

r r-^j 

^H-2e) = P^,„(r), 

2, . 2iuj 2A/ 

—u! hi + iLuhn. ho H Trhi + 

r 



96z7rawM 



^/lo + fh'i + ^^hi = Pu,eni{r) ■ 



(A2) 

(A3) 
(A4) 



where A ~ {t + 2){£ ~ l)/2 and the source terms arc the 
Fourier transforms of 



a -2)1 



(A5) 



Pit, r)^ l&^r^'.^^^^ J T^^WXb dn . (A6) 

The three equations above are not independent, due to 
the Bianchi identity 



1. Perturbation equations with a general source 



In our perturbative approach, the spacetime metric is 
-ti' = + h^n where g'l^} is the Schwarzschild metric 



(EH) 



and hfj^^ is the metric perturbation. We decom- 
pose h^i, in tensor spherical harmonics choosing the RW 
gauge, as in Ref. [2^ where the perturbation equations 
without source were derived. 



There are two families of perturbations, the odd (or 
axial) parity perturbations (in our gauge described by 
the functions /ig™(aj,r) and ft,^'"(w,r)), and the even 
(or polar) parity perturbations (in our gauge described 
by the functions H^"'{uj,r), Hi"^{oj,r), Hl'^{uj,r) and 
K^"^{uj, r)). Moreover, we decompose the scalar field as 



(Al) 

In the rest of the Appendix, we leave implicit the uj de- 
pendence and the i and m indices, a sum over which is 
assumed. The relevant first order Einstein equations for 



dtP' + drP'' + -P'' - ^P - . 



2_A 

r2 



Jtr -r Ur-T -I -T (A7) 

Defining the RW function as (notice the sign difference 
with respect to Ref. (22| ) 



{oj,r) 



(A8) 



from Eqs. (jA2p - (jA4p we obtain the equation for axial 
parity perturbations (|2.9p . with source term given by 



(A9) 



SRw{t,r) = - 
r 



Ui-'-f-)p-fd,.p+p^- 



(AlO) 



The perturbation equation for the scalar field can be com- 
puted replacing Eq. (|Aip into Eq. p.4p and linearizing it. 
Using Eq. (|A3|) . we obtain Eq. (|2.10p with source term 
given by 



Ssit,r)=-if{r) 



6^(^+l)aAf 



pr 



(All) 
(A12) 



We note that this source term is entirely due to the CS 
coupling. The polar sector, instead, is unaffected by 
the CS coupling. Hence, the equation for polar pertur- 
bations (|2.1ip is the same as in general relativity with 
source term given by (32j : 
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SziuJ,r) ^ JdtSz{t,r)e 
Sz{t,r) 



1 



(A + 1)A 
P 



or or 



+ r(A - f)N"' + rfN^ 



A(A - l)r^ + (4A - 9)i\/r + 15M^ 



A r 



(A13) 



(A14) 



where A = A + 3A/ / r and we have defined 



A^^ = 87rr- 



AT" = 



W+T) 



327rr4 



l)£(^+l)(^ + 2) 
I 



(A15) 
(A16) 



2. Source term for point particle on geodesies 

Here we work out the exphcit source terms for a point 
particle moving on geodesies around a Schwarzschild BH 
(see e.g. Ref. [32[)- Let us consider a point- hke particle 
on a timelikc geodesic with coordinates (t) . The stress- 
energy tensor (in the time domain) is 



dr 



i^'u" (a;" - z^) 



(A17) 



where /x is the mass of the particle, S'^{x°' — z^) is the 
four-dimensional Dirac delta, r is the proper time and 
u'^{t) = z^{t) is the four- velocity. The integral above 
can be explicitly computed: 

We introduce the semi-latus rectum p and the eccentric- 
ity e, as orbital parameters. They are defined so that 
the periastron and apastron are at r = pAf/(l -I- e) and 
r — pM / (\ — e), respectively. In terms of these param- 
eters, the energy and angular momentum per unit mass 
of a point particle are 



(p-2-2e)(p-2 + 2e) 



pip-S- e^) 



p — 3 — 
and the four velocity reads 

i E 



(A18) 



(A19) 



where = f [1 + L^/' 

For point-like particles on geodesies the time- 
dependent source terms (|A10[) and (|A14p can be com- 



puted explicitly. They read [32| 

SRw{t, r) = G^wit, r)S[r - rp{t)] + 
+FRw{t,r)S'[r-rp{t)], 
Sz{t,r)^Gz{t,r)S[r-rpit)] + 
+Fzit,r)S'[r-rp{t)], 



(A20) 
(A21) 



with 

GRwit,r) 



A + B 



r / 



J' 



FRw{t,r) = -A—, 

Gz{t,r)^aY*{t) + bY^; {t) + c u;^ (t) + d v;^ it) , 

(A22) 



Fz{t,r) 
where 
A 



^ Y*{t), 



\ + l K E 



167r 



{i^2)\ E '^'^ 



B 



Stt P I 6M ~ 



E 



MI 



A + 1 rA2 

E L r 
167r L p ^ 

11 

\ + lEr^K ' 

2)\P p 



A + 3 

8^ f 
X + 1 E r^A' 



7M 

r 

2 f3 



d = -327r 



{£ + 2)1 E r3 



(A23) 



In the above expressions, Y{t), X^(i), UAsit), VAsit), 
WAsit) denote the scalar, vector and tensor spheri- 
cal harmonics evaluated at the angular position of the 
particle ipp{t); thus, Y{t) is shorthand notation for 
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Y^"^{'K/2,ipp{t)), and the same holds for the other har- 
monics. 

Since the orbital motion takes place in the equatorial 
plane, each spherical harmonic function is evaluated at 
9p = 7r/2. An important consequence of this is that the 
source term for the Zerilli function vanishes when i+m is 
odd, while the source term for the RW function vanishes 
when £ + m is even. 

Finally, the source term for the scalar equation (|A12[) 

is 



Ssit,r) = Gs{t,r)S[r~rp{t)], 



(A24) 



where 
Gs{t,r) 



2 967raM L \J & ~ V"' 



X*(t). (A25) 



Notice that Ss, at variance with Srw and Sz, does not 
contain derivatives of the Dirac delta. Note also that if 
the orbit is circular, E ~ V, thus 5*^ = 0. This can 
be traced back to Eq. (jA12[) . and to the fact that ^ 
T^^ ~ m'' = for circular orbits. 



3. Source describing a particle in circular orbit 

All quantities in the previous section are considered in 
the frequency domain: they depend on r and w, although 
the dependence on lo has often been left implicit. To com- 
pute the sources S'z(w, r) and Sb.w{'^^ t), one should first 
consider the time domain sources Szit, r), Sawit, r), and 
then compute their Fourier transforms. This operation 
proceeds straightforwardly in the case of a circular orbit 
r = f. Indeed, we have 

rp{t)=f, (j>p(t) = cuKt , 

where the Keplerian frequency reads 



The dependence on t comes from the tensor spherical 
harmonics. For instance Eq. (|A22p now reads (hereafter 
we write explicitly the indexes £, m) 



GTit.r)=aY*'"\ep{t),Mt)) + 



(A29) 



with 



r*^™(0p(t),0pW)=y*^™(-,o)e 



i'm<pp{t) 



and similarly for the other terms in Eq. (|A22[) . Hence, we 
can write G'^'(t, r) as a quantity which does not depend 
on t (conventionally we will indicate it with a hat) times 



-imuJKt. 



G'|"(t,r) = G|"(r)e-'™"^^* 



where 



Gfir) =aY*'"' {-,0 



(A30) 



(A31) 



(A26) 



The same can be done with all other quantities. Then, 
the Fourier transform gives: 

G|"(a;,r)=— / d<G|™(r)e-™'^^*e''^* 

= G|™(r)(5(L.-mc^K), 
and the same holds for the other quantities. Therefore, 
S'|j",V(w, r) =S{u}~ mujK) x 

(A32) 

S'|™(w, r) = S{uj - rriLOK) x 

G'|"(r)(5(r-f)-f F|™(r)y(r-f)" . 

(A33) 



Furthermore the geodesies energy, angular momentum 
and four-velocity respectively read 



E 



L = 



2M 




(A27) 



Using the definitions above, the source terms (|A20[) . 
(|MT|) . reduce to 

SRwit, r) = Gw(i, r)S{r - f) + Fnw{t, r)S'{r - f) 
Sz{t, r) = Gz{t, r)S{r - f) + Fz{t, r)5' [r - f) 
Ss{t,r)^{). (A28) 



Appendix B: Perturbative Green's function 
approach 

In this Appendix we develop a perturbative Green's 
function approach to solve Eqs. (|2.9p - (|2.1ip . This ap- 
proach is valid in the small coupling limit. We shall com- 
pare it with the general method (valid for any coupling) 
adopted in the main text. The perturbative approach can 
be useful for possible analytical calculations and it is also 
important as independent check for numerical results. 

As long as the dimensionless CS coupling C (|2.19p is 
small, Eqs. ()2.9p - ()2.11|) can be solved by a perturbative 
scheme. Indeed, as we shall see, the Regge- Wheeler func- 
tion and the scalar field are consistent with the ansatz 



(Bl) 
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The function Q^^^ is the solution of the standard Regge- 
Whceler equation in general relativity 



Q{r) = Sfiwir) 



(B2) 



and can be solved by the Green's function approach, as 
in Section lill Al by considering the particle source only: 



g(°)(r) 



Q+jr) 
Wrw 



■I dr^.Q_Sjiw + -rrr / dr^Q+S_ 

J~oo yyRwJr 



RW 



(B3) 



(note that the Green's function expression is linear in the 
source). The function G^*^-* is the lowest order contribu- 
tion (in the CS coupling) to the scalar field 9. Its equa- 
tion can be found replacing the ansatz (jBip in Eq. (|3.6p : 
at lowest order in C, it gives 



dr 



— +uj'-Vs 



e 



(0) 



Ss-f 



{£ + 2)l6iMa 
{£-2)1 wr5/3 ^ 



)(0) 



(B4) 

where Ss ^ a/ /3 is given in Eq. (jA12p and Vs is Eq. (j2.14p 
at order a. Defining an effective source term 



(B5) 



we can write the solution of (IB4I) as 



e+(r) 



r+oo 

+e_(r) / dr,e+Ss 



(B6) 



Finally, if we denote the scalar field source in the Regge- 
Wheeler equation (at lowest order in C) as 

SRwir) = ae^^(r), (B7) 

we get that the main correction to the Regge- Wheeler 
function is 



CQ 



(1)- 



1 



Wrw 



Q+{r) 

+ 00 



dr^,Q^Siiw 



dr^:Q+SRw 



(B8) 



Note that Ss ^ therefore G*-"-* ^ oi/ 13, and then 

Srw ~ ae(°) ~ C, consistently with Eq. (|Bg)) . 

At infinity and at the horizon, where the fluxes are 
computed, we get 

g(n ^ ±oo)= (Q(o) + CQ^'^)(n ^ ±oo) 

dr^Q^{SRw + Srw) ■ 



RW 



(B9) 



Even for circular orbits, Srw ~ B^°^ is defined through- 
out the entire spacetime, thus the evaluation of the in- 
tegrals above is more involved than in general relativity, 
in which we simply have Srw ^ ^if — ^) for a circular 
orbit. This makes the perturbative approach much more 
time-consuming than the non-perturbative approach. 

It is straightforward to extract the quantities relevant 
to compute the fluxes, Eqs. (|4.2p and (|4.4p . Indeed, by 
applying the same procedure explained in the main text, 
we can deflne Q^^^™(a;,r) = Q^^^™(r)(5(aj — mw/^)e^*"''* 



with 



:(r) 



Wrw 



Q±G 



RW 



Q±F 



RW 



f 



f 



for r ^ r respectively. Then, from Eqs. (jB5])-(jB8j) one can 
compute the quantities 6'j?''^™(r), 5'^^(r), Q£^^'"(r), 
and finally the complete Regge- Wheeler function Qj!" = 



Q{0)£m 



f C2± sit infinity, which has the form 
Q^™(c^) = Q'^6{uj - mujK)e^''^'-' . (BIO) 



In the same way, the scalar field perturbation 

e 



(w) = 9j_ (a;,r* ±oo) has the form 



(0)£ri 



(^) 



'S{uj — mLdK)e 



(Bll) 



Notice that 6^"^ ~ a//?, so that the corresponding energy 
flux is 0{(), consistently with our expansion (jBSp . The 
contribution of (Q^^^ to the energy flux is ©(C^) and it 
is neglected in this approximation. 



a. Comparison between the perturbative and the 
non-perturbative approach 

Comparing the perturbative method with the general 
approach discussed in the main text is important for two 
reasons. First, it gives an independent check of both 
methods. Second, it allows to assert the validity region 
of the perturbative approach. In Fig. [S] we compare the 
corrections to the total energy flux and to the number of 
cycles obtained by the two methods. 

When C <C 1, the methods agree very well, confirm- 
ing each other. For example, when p ~ 6, the emit- 
ted fluxes computed with the two methods agree within 
0.1% or better for C < 0.3, whereas they differ by 
~ (1%, 5%, 10%, 50%, 100%) when C ~ (0.5,1,2,10,20), 
respectively. Remarkably, the perturbative approach is 
valid within a few percent up to C ^ 1 whereas, for larger 
values of C, it gives an overestimated flux. 



b. On the dominance of fluxes at the horizon in DCS theory 

One of the curious results borne out of this study, 
with important consequences for the main observational 
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V m^/M^ 

FIG. 5. (Color online) Comparison between the non-perturbative method described in the main text and the perturbative 
method described in Appendix [B] Left: Corrections to the total flux emitted. Right: Corrections to the number of cycles 
for 1712 ~ IAMq. In both panels, from below to top: = 0.01,0.1,0.5,1,2,5,10,20. For C 1 the perturbative method 
overestimates the total flux, resulting in a larger (in absolute value) correction to SM. 
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FIG. 6. (Color online) Homogeneous solutions of the Regge- 
Wheeler equation, \Q±\. 



prospects, concerns the fluxes at the horizon: while in 
general relativity fluxes at the horizon are orders of mag- 
nitude smaller than at infinity, this does not occur in 
DCS gravity; as summarized in Eqs. (|4.8p the corrections 
imparted by the DCS coupling affect more strongly the 
flux at the horizon rather than at infinity. This can be 
understood from our perturbative analysis in the follow- 
ing way. From the above study, the axial Regge- Wheeler 
function in general relativity behaves as 



)(0) 



(r* — > ±oo) 



Wrw 



+ CXD 



di\Q^SRw 



(B12) 



We assume for simplicity circular orbits and that Srw ^ 
6(r ~ f), i.e. we neglect the S'{r — f) contribution; how- 
ever, this discussion can be easily generalized to non- 
circular orbits. The fluxes at the horizon and at infinity 
then depend on the value of the homogeneous solutions 
at r = f. From Fig. |6] where we plot these homogeneous 
solutions, we see that, for r > f > 6M, 



\Q- 



\'>\Q-\' 



(B13) 



which explains why the flux at the horizon is smaller than 
that at inflnity. On the other hand, in this perturba- 
tive approach, the DCS correction to the Regge- Wheeler 
function behaves as 



CQ''^\r, ^ ±oo) 



1 



Wr 



w 



-t-oo 



dr^,QzrSjiw 



B14) 

Since SRw{r) Q^^^r^, the largest contribution comes 
from the near-horizon region 2M < r < 6M, where (see 
FigiD 



IQ+I < IQ- 



(B15) 



therefore the largest DCS correction to the flux is at the 
horizon. 

Analogously, the DCS scalar field behaves as 



e("'(r, ^ ±oo) 



1 



Wrw 



+ 00 



dr^Q^Ssir) 



(B16) 

The homogeneous solutions of the scalar equation, Q±, 
are similar to the Regge- Wheeler homogeneous solutions 
Q±; in particular, close to the horizon one has |e+r < 
|e_|^ and Ss{r) - Q(°V^^ therefore the DCS scalar 
flux is larger at the horizon than at infinity. 
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Appendix C: Comparison of different ways to 
estimate the frequency one year prior to merger 

In our calculation of the evolution of orbital frequency 
with time, we have adopted an adiabatic approxima- 
tion, where the fluxes are determined numerically from 
those of a particle in geodesic motion. In the literature 
one sometimes finds other alternative calculations, with 
which we now compare our (more accurate) results. 

One of the most common alternatives for estimating / 
in Eqs. (j4.9p and (|4.12p consists in taking a PN approxi- 
mation which, at 2PN level, reads l43l : 



OTT 



/34103 13661 



59 



V 18144 2016 
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(CI) 



(see e.g. Refs. [50, |5l| for the 3.5PN formula). In the 
equation above, M = mi + TO2 is the total mass of the 
two-body system and rj = mim2/M^ = ji/M . Neglecting 
2PN terms in Eq. JCl]), we can solve Eq. for /ly^ 

analytically: 



/: 



lyr 



5^/8 
2y/2'KM 



405 + 16?7 



Jobs 



-3/8 



(C2) 



For example, from the formula above /lyr ~ 0.00274 Hz 
for Tobs = lyr, mi = lO^Mg and TO2 = IOMq, while, 
including 2PN corrections, for the same parameters we 
obtain /urr 0.00252 Hz, and including 3.5PN cor- 
rections [50, [EH we get /lyr ^ 0.00242 Hz. Our own 
geodesic-based approach plus the numerical fluxes yields 
/lyr ~ 0.00223 Hz. Notice also that, if we insert the 
quadrupole formula for the energy flux of a particle in a 
circular orbit, 



E = E 



N 



32 ^^A/S 



(C3) 



in Eq. (|4.16p and expand for r/M ^ 1, then we recover 
the PN formula (ICll) at first order. 



The errors introduced by the PN approximation may 
be evaluated with a hybrid approach as follows. We com- 
pute the modification SAf using the numerical flux E^cs 
in Eq. (|4.19|) but, in order to estimate the lower boundary 
of the integral (|4.12p . we use four different prescriptions: 
(i) / is computed numerically within our geodesic-based 
approach, as explained in the main text; (ii) / is given 
by its 3.5PN expansion (see Eq. (32) in Ref. [5l|); (in) / 
is given by its 2PN expansion (jCl[) : (iv) / is given by the 
truncation of Eq. (jCl[) at Newtonian order. These dif- 
ferent prescriptions affect the value of /in and, in turn, 
Eq. (I4l9l) . 

This is shown in the left panel of Fig. [T] For large 
central masses, i.e. low frequencies, all the different pre- 
scriptions yield basically the same result to a good accu- 
racy. For smaller central masses, we expect relativistic 
effects to become important: the system enters the LISA 
band when the small mass is already close to the ISCO, 
and where the PN expansion is less accurate. Indeed, in 
this regime the PN formula can give factors of order ~ 2 
difference with respect to the more accurate prescription 
adopted in the main text. 

As shown in the right panel of Fig. [3 similar devia- 
tions are observed also in general relativity. Indeed, if 
we compute Af from Eq. (j4.9|) using different prescrip- 
tions for /in, we find that the number of cycles computed 
using the PN formula may be overestimated by ^ 20%. 

As discussed in Section lIVBi the method based on 
the adiabatic approximation is better suited to deal with 
EMRIs. We note here that this is true not only in al- 
ternative theories like DCS gravity, but also in general 
relativity. 



Appendix D: Expressions for the coefficients Cj_ 



(i) 



and D 



(0 



The coefficients in Eqs. p.l3p and p.l4p depend solely 
on the solutions of the homogeneous system associated to 
Eq. (|3.7p and their explicit forms read 



C 



(1) 



A" 



0(i)ga)e'(2) 



0a)Q(2)QKi) _ 0(i)Q(2)0/a) _ e(i)Q(i)e,(2) 



0(2)^(1)0,(1) 



c_ 



(1) 



c 



(2) 



A-1 
A-i 



0(l)g(2)0,(2) _ 0(^2)^(2)0,(1) _ 0(2)q(1)0,(^2) ^ ^(2) (2) ^ 0(2)^(2)0,(1) _ Qi2) ^,(2) 

-e^gWe'f + e^^gL^'e'^!' + eL'^gi'^e'f - g^e'W - eL^'g^^e'^ + e^gf el^^l (pi) 
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FIG. 7. (Color online) Left: Same as right panel of Fig. (2] but with different prescriptions to compute the frequency one year 
prior to merger. The first prescription ("GR") is adopted in the main text and takes circular geodesic motion and the energy 
fluxes numerically computed in general relativity for a particle in circular orbit. The second ("3.5PN") and third ("2PN") 
prescriptions take respectively the 3.5PN formula [sol . and the 2PN formula, Eq. (|C1[I for the evolution of the frequency, 
while the third ("Newtonian") is its truncation at lowest order. Right: relative difference between the total number of cycles 
in different approximations with respect to our geodesic-based approach for ^ — (general relativity). At low frequencies all 
prescriptions yield identical results, while at higher frequencies (smaller central masses), they can differ by a factor of order 
two. 



and 



(2) 



D 



where 



A- 



-ef g'L'^gW + eV'^g'L'^gf + g'V'^gL'^ - eL'^g'V'^gf - e^'^g'f gL^' + eL'^g'f g« 
e?)g'L^^g(_^) - e^^^^g'L^^g?) - e^^g'V^'gL^) + eL^^g'V^'g?' + e^^^^g'f gL^' - eL^^g'f g« 



,(D2) 



A = e 



/(2) 



/(I) 



g 
g'f 



ea)QKi)Q(2) _ eL^^g'WgV') - e^gl'^gL^^ + eL^^g'L'^gV'^' 

0(2)^(1)0,(2) _ 0(^2)^(1)0/(2) _ 0a)Q(2)0,(2) ^ 0(^2)^(2)0,(1) ^ 0(1) q(^2) 0,(2) _ 0(2) q(2) 0,(1) 

-eL^'gL^'e'W + eV^'gL^^e'L^' + e^gL^'e'^ - e(.^'gL')e'L^' - e^gV^^e'L'^ + eL^'g^'^e'^ 
ef g'WgL^'e'W + g'L'^gV''e'L'^ + eL^'g'Wg^^e'^ _ e^'^g'L'^gf e'L'^ + ef g'L'^gL^^e'V'^ 



(D3) 



r 



is a generalized Wronskian, it is constant by virtue of out of the integrals in Eqs. p.l3p and p.l4p . In the 
the homogeneous system, and therefore can be factored equations above a prime denotes derivative vifith respect 

to the tortoise coordinate r,. 
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